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Abstract 

We calculate the branching ratios of pure annihilation type decays B° — > D~K^ + and 
B s — > Dai using the perturbative QCD approach based on kx factorization. The branching 
ratios are predicted to be (60.6±^j ±i 3 4 + - 2 2 i) x for B° D~ K* + , (l-lioiloita^xlO" 6 
for B s Dal and ( 2 - 3 ias ^aliai) x 10 ~ 6 for B * -> D~ a t ■ Tlie Y are lar S e enough to be 
measured in the ongoing experiment. Due to the shortage of contributions from penguin 
operators, there are no direct CP asymmetries for these decays in the standard model. We 
also derive simple relations among these decay channels to reduce theoretical uncertainties 
for the experiments to test the accuracy of theory and search of new physics signal. 
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I. INTRODUCTION 



Two body hadronic B decays have been a hot topic for many years, since it involves 
the perturbative QCD calculation and factorization study. It is also important for the 
test of standard model, the CKM angle measurements and the search of new physics 
phenomena. For many years, people do calculations based on the naive factorization 
— on, Me, proved* the soft-coHinear effective theorvfl . However, there is one 
kind of diagrams, the so called annihilation type diagrams, which was argued to be 
helicity suppressed since no one knows how to calculate. In the well developed collinear 
factorization, there is endpoint singularity in the calculation of these diagrams. In fact, 
this kind of diagrams are essential for the strong phase of direct CP asymmetry in the 
B — > K + tt~ decays [2], which is proved to be important. 

Furthermore, there is one kind of B decays, which contains only annihilation type di- 
agram contributions. One of the examples is the B° — > D S K + decay, which is predicted 
in ref. 3-5] and measured by the B factories later jfj. Recently, the CDF collaboration 
measured the first pure annihilation type decays in the B s sector i.e. B s — > 7r + 7r _ de- 
cay, which exactly confirms the perturbative QCD prediction for this decay [7, 8|. It is 
worth of mentioning that the perturbative QCD (PQCD) approach is almost the only 
method can do the quantitative calculations of the annihilation type diagrams 7|, |8|. 

In this paper, we shall study the pure annihilation type charmed decays B° — > 
D~K2 + and B s — > Da<i in the PQCD approach, which is based on the kx factorization 
[q|,[io|. These decays are predicted to have a large branching ratio as 1CT 6 to 1CT 5 , which 
are measurable in the near future experiments. In the annihilation type diagrams, both 
of the light quark and the heavy anti b quark in B meson annihilate into another quark 
anti-quark pair through the four quark operators, while another light quark pair in the 
final state mesons are produce by a gluon attaching to the four quark operator. Since the 
light quark in the final states are collinear, the gluon connection them must be hard. So 
the hard part of the PQCD approach contains six quarks rather than four quarks. This 
is called six-quark effective theory or six-quark operator. In this approach, the quarks' 
intrinsic transverse momenta are kept to avoid the endpoint divergence. Because of the 
additional energy scale introduced by the transverse momentum, double logarithms will 
appear in the QCD radiative corrections. We resum these double logarithms to give 
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a Sudakov factor, which effectively suppresses the end-point region contribution. This 
makes the PQCD approach more reliable and consistent. 

This paper is organized as following. In Sec. II, we present the formalism and perform 
the perturbative calculations for considered decay channels with the PQCD approach. 
The numerical results and phenomenological analysis are given in Sec. III. Finally, Sec. IV 
contains a short summary. 



II. FORMALISM AND PERTURBATIVE CALCULATION 

The B° — > D~K^ ,B S — > D°a\ and B s — > D~a\ decays are pure annihilation type 
rare decays. At the quark level, these decays are described by the effective Hamiltonian 



H eff 
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H eff = ^ v; b V uD [dMOiQi) + C 2 (/i)0 2 (/i)] , (1) 

where V c b and V u d are CKM matrix elements, U D" denotes the light down quark d or 
s, and Ci )2 (yu) are Wilson coefficients at the renormalization scale fi. Oi i2 (/i) are the 
four quark operators. 

0\ = (b a Cp)y_ A (upD a )y_ A , 2 = {b a C a )v-A{UpDp)y_ A . (2) 

where a and are the color indices, (baC^v-A = ^7^(1 — 7 5 ) c /3- Conventionally, we 
define the combined Wilson coefficients as 

a t = C 2 + Ci/3, a 2 = C 1 + C 2 /3. (3) 

In the hadronic matrix element calculation, we factorize the decay amplitude into 



soft(<I>), hard(H), and harder (C) dynamics characterized by different scales [12l. Il3|. 



A ~ J dx\dx 2 dx-ib\db]b 2 db 2 b-idb-i 

xTr [C(t)^ B (x u b 1 )^ M2 (x 2 ,b 2 )^ M3 (x ;i ,b 3 )H(x h b u t)S t (x l )e- s{t) ] . (4) 

where bi is the conjugate variable of quark's transverse momentum kiT, Xi is the momen- 
tum fractions of valence quarks, and t is the largest energy scale in function H(xi,bi,t) 
which is the hard part. C(t) are the Wilson coefficients with resummation of the large 
logarithms \\i{m w /t) produced by the QCD corrections of four quark operators. S t (xi) 




is the jet function, which is obtained by the threshold resummation and smears the 
end-point singularities on Xi 14]. The last term, e~ s<yt \ is the Sudakov form factor, 



from resummation of double logarithms, which suppresses the soft dynamics effectively 



and the long distance contributions in the large b region |15l. |16|. Thus it makes the per- 
turbative calculation of the hard part H applicable at intermediate scale, i.e., m B scale. 
The $j, meson wave functions, are nonperturbative input parameters but universal for 
all decay modes. 

The lowest order Feynman diagrams of the considered decays are shown in Fig.l. 
The amplitude from factorizable diagrams (a) and (b) in Fig.l is 



[2 f 1 f 1/A 

A af = 8y -C F f B ^m B J dx 2 dx 3 J b 2 db 2 b 3 db 3 (f) D (x 2 ,b 2 ) 

x {[-<f) T {x 3 )x2 + 2r D r T (f) s T (x 3 ){x2 + 1)] 
■h af ((l - x 3 ), x 2 (l - r 2 D ), b 2 , b 3 )E af (t e ) 

- [0t(x 3 )(x 3 - 1) - r D r T (4(x 3 )(l - 2x 3 ) + <p s T {x 3 ) (2x 3 - 3))] 

■ h af (x 2 , (1 - x 3 )(l - rf,), 63, b 2 )E af (t f )} . (5) 

In this function, Cp = 4/3 is the group factor of SU(3) C , and r D ^ = m D ^/m B . The 
hard scale t e j and the functions E a f and h a f are given by 

t e = max{^/x 2 (l - r 2 D )m B , l/b 2 , l/6 3 }, 

t f = max{^/ (1 - x 3 )(l - r 2 D )m B , l/b 2 , l/h}, (6) 



E af (t) = a s (t) ■ exp[-S T {t) - S D (t)}, (7) 

A 



h af (x 2 ,x 3 ,b 2 ,b 3 ) 



<"S2o-(l) 



9{b 2 - b 3 )H^> {y/x3~m B b 2 ) Jo {\fx2,m B b^ + 
0(b 3 - b 2 )H { 1) {y/x^m B b 3 ) J {y/x^m B b 2 ) ■ S t (x 3 ). 



The amplitude for nonfactorizable diagrams (c) and (d) in Fig.l is 

32 f 1 f 1 ^ 

Manf = -irC F 7rm% / dx x dx 2 dx 3 / bidbib 2 db 2 <f) B (xi,bi)<f) D (x 2 ,b 2 ) 
<J Jo Jo 

x { [<p T {x 3 )x 2 + r D r T ((j) s T (x 3 ){x 3 - x 2 - 3) + 4> t T {x 3 )(x 2 + x 3 - 1))] 
■h an fi(x 1 ,x 2 , x 3 , 61, b 2 )E anf (t g ) 

+ [<fo(x 3 )(x 3 - 1) + r D r T (^(x 3 )(x 2 - x 3 + 1) + </4,(x 3 )(x 2 + x 3 - 1))] 
• h anf2 (x 1 , x 2 , x 3 , bi, b 2 )E anf (t h )} . (9) 



t 9 = max{yz 2 (l - x 3 )(l - rjj)m B , yjl - (1 - (1 - rr 3 )(l - rf>))(l - xi - x 2 )m B , 
1/61,1/62}, 



t h = max{y x 2 (l - x 3 )(l - r 2 D )m B , yj (1 - ar 3 )(l - rf,)|<Pi - x 2 |m B , 

1/61,1/62}, (10) 



E an f = a s (t) ■ exp[-S B (t) - S T (t) - S D (t)] \ b2 =b 3 , 
{x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) 



with j — 1,2. 
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x -b 2 )H^ {Fm B b x )J {Fm B b 2 ) 

+6(b 2 - b x )H^ (Fm B b 2 ) J {Fm B b x ) 
(1) 



ITT TT\ 



x < * u V V 1 J 
if (Fjirisbi) 



\F?\m B b x ), F?<0, 



if >0, 



F 2 = x 2 (l-x 3 )(l-r 2 D ), 

F 2 = 1 - (1 - (1 - x 3 )(l - rf,))(l - X! - x 2 ) ; 

F| = (l-ar 3 )(l-7i,)(xi-x 2 ). 

The expressions of S B (t), Sx(t), Sr>(t) and £4 can be found in ref. 
wave functions of initial and final states can be found in ref. 14 
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(13) 
The 



With the functions obtained in the above, the amplitudes of these pure annihilation 
type decay channels can be given by 

A(B° -+ D;K* 2 + ) = ^Zv; b V ud [a 2 A af + C 2 M mf ], (14) 



A(B° S D°a° 2 ) = ^-Lv; b V us [a 2 A af + C 2 M anf ], (15) 



A(B° S -+ D~<4) = ^V; b V us [a 2 A af + C 2 M anf }. (16) 



III. NUMERICAL RESULTS AND DISCUSSIONS 

For numerical analysis, we use the following input parameters: 

f B/Bs = 0.21/0.23 GeF, f D/Ds = 0.205/0.241 GeV, = 118(77) MeV, 
/ip = 102(117) MeV, M D/Ds = 1.869/1.968 GeV, M B/Bs = 5.279/5.366 GeV, 
\V cb \ = 0.0415 ± 0.0011, \V ud \ = 0.9742 ± 0.0002, \V US \ = 0.2257 ± 0.0012, 
A f = c 4 D = 0.25 GeV (17) 

After numerical calculation, the branching ratios of these decays are: 

Br(B° -+ d;k;+) = (eo.el^j^ 2 ,) x io" 6 , 

Br(B s -+ D°a° 2 ) = (l.llg^^i) * IO" 6 , 

Br(B s D-at) = (2.3^j^l°;l) x 10- 6 . (18) 

The branching ratio obtained from the analytic formulas may be sensitive to many 
parameters especially those in the meson wave function. For the theoretical uncer- 
tainties in our calculations, we estimated three kinds of them: The first errors in our 
calculations are caused by the hadronic parameters, such as the decay constants and 
the shape parameters in wave functions of charmed meson and the B^ 8 ) meson, and the 
decay constants of tensor mesons. The second errors are estimated from the unknown 
next-to- leading order QCD corrections with respect to a s and nonperturbative power 
corrections with respect to scales in Sudakov exponents, characterized by the choice of 
the Aqcd = (0.25 ± 0.05) GeV and the variations of the factorization scales defined 



in eq]6] and eqJTOl The third error is from the uncertainties of the CKM matrix ele- 
ments. It is easy to see that the most important theoretical uncertainty is caused by 
the non-perturbative hadronic parameters, which are universal and can be improved by 
experiments. 

These pure annihilation type decays considered in this work are dominant by W ex- 
change diagram. All these decays do not have contributions from the penguin operators. 
Since the direct CP asymmetry is caused by the interference between the contributions 
of tree operators and that of penguin operators, it does not appear in these modes. 
Although the annihilation type diagrams are power suppressed in PQCD approach, the 
branching ratio of these considered Cabibbo-Kobayashi-Maskawa-favored decays are 
sizable and large enough to be measured in experiment. Through the study of these 
pure annihilation type decay modes, we can understand the annihilation mechanism in 
B physics well. 

It is easy to find that there are large theoretical uncertainties in any of the individual 
decay mode calculations. However, we can reduce the uncertainties by ratios of decay 
channels. For example, simple relations among these decay channels are derived from 
eq.tdUHSD 

BrjB -> D-K;+) 2f D yl A 60.6 

— r^j rsj nil r^j 

Br(B s D«al) f*V* 1.1 ' 

Br(B s -> D°a° 2 ) ^ 1 ^ LI 

Br{B s D-a+) ~ 2 ~ 2.3' 1 ' 

It is obvious that any significant deviation from the above relations will be a signal of 
new physics. 

IV. SUMMARY 

We calculate the branching ratios of three pure annihilation type decays in the 
perturbative QCD approach. The predicted branching ratios are Br(B° — > D~K 2 ) ~ 
6 x 10~ 5 , Br(B s ->■ D°a° 2 ) ~ 1 x KT 6 and Br(B s ->■ D'at) ~ 2 x KT 6 . They are 
sizable and large enough to be measured in forthcoming experiment. The study about 
the pure annihilation type decays can help us understand the annihilation mechanism 
in B physics. There are no direct CP asymmetries, because these decays have no 
contributions from penguin operators in the standard model. 
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